On the equilibrium of self-gravitating neutrons, protons and electrons in /3-equilibrium 
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We have recently proved the impossibility of imposing the condition of local charge neutrality in a 
self-gravitating system of degenerate neutrons, protons and electrons in /J-equilibrium. The coupled 
system of the general relativistic Thomas-Fermi equations and the Einstein-Maxwell equations have 
been shown to supersede the traditional Tolman-Oppenheimer-Volkoff equations. Here we present 
the Newtonian limit of the new equilibrium equations. We also extend the treatment to the case 
of finite temperatures and finally we give the explicit demonstration of the constancy of the Klein 
potentials in the case of finite temperatures generalizing the condition of constancy of the general 
relativistic Fermi energies in the case of zero temperatures. 



I. INTRODUCTION 

It is well known that the classic work of Oppenheimer 
and VolkofF T| addresses the problem of neutron star 
equilibrium configurations composed only of neutrons. 
For the more general case when protons and electrons 
are also considered, in nearly all of the scientific litera- 
ture on neutron stars it is assumed that the condition of 
local charge neutrality applies identically to all points of 
the equilibrium configuration (see e.g. (4I). Consequently, 
the corresponding solutions of the Einstein equations for 
a non-rotating neutron star, following the work of Tolman 
3 and of Oppenheimer and Volkoff [i|, have been sys- 
tematically applied. However, the necessity of processes 
leading to electrodynamical phenomena during the grav- 
itational collapse to a black hole [J] suggests a critical 
reexamination of the current treatment of neutron stars. 
In this regard in a set of recent articles (see e.g. [a, Q) 
we have developed the first steps toward a new consistent 
treatment for the description of neutron stars overcoming 
the traditional Tolman-Oppenheimer-Volkoff equations. 

We have recently proved the impossibility of impos- 
ing the condition of local charge neutrality in a self- 
gravitating system of degenerate neutrons, protons and 
electrons in /3-equilibrium [y] . It has been there discussed 
that the traditional approach for the description of neu- 
tron stars adopting the condition of local charge neutral- 
ity is inconsistent with the Einstein-Maxwell equations 
and with micro-physical conditions of equilibrium within 
the quantum statistics. We emphasize the basic role of 
the constancy of the general relativistic Fermi energy 
of each species. Consequently, the traditional Tolman- 
Oppenheimer-Volkoff system of equations has been su- 
perseded by the coupled system of the general relativis- 
tic Thomas-Fermi equations and the Einstein-Maxwell 
equations. We have been guided in this new approach 
by the generalization of the Feynman-Metropolis- Teller 
treatment of compressed atoms to the relativistic regimes 
0. There, using the existence of scaling laws 043; ^^^ 



results were extended from heavy nuclei to the case of 
globally neutral nuclear matter cores of stellar dimen- 
sions with mass numbers A « {mpianck/fnn)^ ~ lO^'' or 

Afeore ~ Mq. 

In the present article we first recall in Sec. |n]the new 
system of general relativistic Einstein-Maxwell- Thomas- 
Fermi governing a self-gravitating system of neutrons, 
protons and electrons in /^-equilibrium. We then address 
three new aspects of the problem: 1) we study the New- 
tonian limit of the equilibrium equations (see Sec. Ill); 2) 
we consider the proper generalization of the above treat- 
ment to the case of finite temperatures (see Sec. IV); 
3) we generalize and prove the constancy of the Klein 
potentials of each species generalizing the condition of 
constancy of the general relativistic Fermi energies en- 
countered in the case of zero temperatures (see Sec.lIVp. 
We finally show that the thermal effects do not affect the 
electrodynamical structure of the equilibrium configura- 
tions and the results are qualitatively and quantitatively 
quite similar to the ones obtained with the degenerate 
approximation (see Sec. lIVI) . 



II. EINSTEIN-MAXWELL-THOMAS-FERMI 
EQUATIONS IN THE DEGENERATE CASE 

Following p] we consider the equilibrium configura- 
tions of a degenerate gas of neutrons, protons and elec- 
trons with total matter energy density and total matter 
pressure 
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where ei{p) — -^/c?p^ + m|c? is the relativistic single par- 
ticle energy and Pf denote the Fermi momentum, related 
to the particle number density Ui by rii — {Pf)^/{37T'^h^). 
Introducing the metric for a spherically symmetric 



non-rotating configuration 

the extension to general relativity of the Thomas-Fermi 
equilibrium condition on the generalized Fermi energies 

Ef = e^/^fM - m,c^ + q,V , (4) 

{i — n,p,e, Qi — 0, e, — e) and the condition of /3- 
equilibrium between neutrons, protons and electrons 
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the full system of equations composed by the Einstein- 
Maxwell- Thomas-Fermi equations can be written as (see 
for details) 
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where e is the fundamental charge, a is the fine struc- 
ture constant, V is the Coulomb potential, /ii = 
dS/dui — y(?\P[Y+rri^l? is the free-chemical poten- 
tial of particle-species, A(r) is the metric function re- 
lated to the mass M{r) and the electric field E{r) — 
— e^'^'^+'^^/^F' (a prime stands for radial derivative) 
through 



l-'-^^Ar'E^r). 
c r c 
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and V ^E[ + eV. 

The condition n^ = Up often adopted in literature is 
not consistent with Eqs. ([5]) and © (see Fig. [T|) there- 
fore we consider the equilibrium configurations fulfilling 
only global charge neutrality (details are given in (6|). 
We solve self-consistently Eq. ^ and (O for the metric, 
Eqs. (|5))- (IT(I1) for the equilibrium of the three degenerate 
fermion species and for the /3-equilibrium. The crucial 
equation relating the proton and the electron distribu- 
tions is then given by the general relativistic Thomas- 
Fermi equation (jlip . The boundary conditions are: for 
Eq. ^ the regularity at the origin: A/(0) = 0, for 



Eqs. ((8])- (fT0|) a given value of the central density, and 
for Eq. (ITT]) the regularity at the origin ne(0) = ^^(0), 
and a second condition at infinity which results in an 
eigenvalue problem determined by imposing the global 
charge neutrality conditions 

ViRe)=E[, T/'(i?e)=0, (13) 

at the radius Re of the electron distribution defined by 

Pf(i?e)=0, (14) 

from which follows 
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Then the eigenvalue problem consists in determining the 
gravitational potential and the Coulomb potential at the 
center of the configuration that satisfy the conditions 
(IT51) - P^ at the boundary. The solution for the den- 
sity, the gravitational potential and electric potential are 
shown in Fig. ([2]) for a configuration with central den- 
sity p(0) = 3.94pnuc- One particular interesting new 
feature is the approach to the boundary of the config- 
uration: three different radii are present corresponding 
to distinct radii at which the individual particle Fermi 
pressures vanish. The radius R^ for the electron compo- 
nent corresponding to Pf{Re) — 0, the radius Rp for the 
proton component corresponding to Pp [Rp) = and the 
radius Rn for the neutron component corresponding to 
P^{Rn) = 0. For the configuration of Fig. [2] we found 
Rn ~ 12.735 km, Rp ~ 12.863 km and R^ - Rp + lO^Ae 
where A^ = H/{mec) denotes the electron Compton wave- 
length. We find that the electron component follows 
closely the proton component up to the radius Rp and 
neutralizes the configuration at Re without having a net 
charge, contrary to the results eg in |ig |. 



III. NEWTONIAN LIMIT 

Despite the fact that the strong gravitational field of 
neutron stars requires a general relativistic treatment, it 
is interesting to explore the Newtonian limit of all the 
above considerations. This can help to elucidate if the 
gravito-electromagnetic effects we have found are of gen- 
eral relativistic nature or to prove their validity in a New- 
tonian regime. 

The Newtonian limit of the equilibrium equations can 
be obtained by the weak-field non-rclativistic limit. We 
expand the gravitational potential at first-order e''" ~ 
1 + $, where the Newtonian gravitational potential is 
<i>(r) — h'{r)/2. In the non-relativistic mechanics limit 
c ^^ oo, the particle chemical potential becomes fii ^■ 
fli + rriiC^, where jli — {Pf)'^/{2mi) denotes the non- 
relativistic free-chemical potential. Applying these con- 
siderations, the electron equilibrium law ^ becomes 
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FIG. 1: Fermi energies for neutrons, protons and electrons 
in units of the pion rest-mass energy for a locally neutral 
configuration with central density p(0) = 3.94pnuc, where 
Pnuc = 2.7 X 10^* g cm~^ denotes the nuclear density. 



which is the classical condition of thermodynamic equi- 
librium in presence of external gravitational an electro- 
static fields (see e.g. [ll[ and [l^]) applied to a gas of 
electrons. 

In the weak-field non-relativistic limit, the Einstein- 
Maxwell equations (|6|- (fTT|) become 
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where p in this case is the rest-mass density 
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The solution of Eqs. dH]), ([T7l)-(l20l) together with the 
/3-equilibrium condition 
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leads to qualitatively similar electrodynamical properties 
as the one obtained in the general relativistic case. In 
Fig. m we show the electric field in the region r < i?„ 
(i?J^°'"* < i?^^) both for the Newtonian as well as for the 
General Relativistic configuration for the given central 
density p(0) = 3.94p„uc- From the quantitative point 
of view, the electric field of the Newtonian configuration 
is larger than the electric field of the general relativistic 
configuration. 



IV. FINITE TEMPERATURE EFFECTS 

The above results have been obtained within the zero 
temperature approximation. It is worth to recall that 
temperatures of the order of ~ 10^ K expected to exist at 
the surface of old neutron stars [ij, [ij] , or temperatures 
of 10^ — 10* K which could, in principle, exist in neutron 
star interiors, should not affect the considerations here 
introduced. For neutron stars, the Fermi temperature 
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where k is the Boltzmann constant, can be as large as 
^ 10^2 K for electrons, ~ 10^^ K for protons and ~ 10^'^ 
K for neutrons for typical central densities of neutron 
stars, which means that neutron star interiors are, at a 
high degree of accuracy, degenerate systems. However, 
the thermal energy associated to temperatures T ^ 10* 
K, EtY, ~ (47r/3)i?3aT4 - 10^* erg, is much larger than 
the Coulomb energy Eq '^ {1/(!>)R^E'^ ^ 10^^ erg, given 
by the internal electric field here considered (see Fig. |4]) . 
Here R is the radius of the configuration and a = a /c 
being a the Stefan-Boltzmann constant. It can be then 
of interest to ask the question if our electrodynamical 
structure will still occur in presence of thermal effects. 

In this more general case, the equation of state given 
by Eqs. (H]) and ^, must be abandoned and replaced by 
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is the Fermi-Dirac fcrmion distribution function which 
gives the particle number density Ui 
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where, for numerical purposes, the free sin gle particle 



energy ei{p) = u{p) - rriiC' 
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rriiC and 

the free particle chemical potential /I^, have been defined 
subtracting the particle rest mass-energy. 

Wc turn now to demonstrate the constancy of the Klein 
potentials throughout the configuration. The equation of 
state ([24|) -([25 |) satisfies the thermodynamic law 
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where s — S/V is the entropy per unit volume and pi = 
dS/drii is the free-chemical potential of the i-specie. At 
zero-temperature T = 0, pi = \J{PfY + mf 
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FIG. 2: Left panel: particle number density approaching the boundary of the configuration in units of the nuclear density 
rinuc — 1-6 X 10"^* cm"^. Right panel: proton and electron Coulomb potential in units of the pion rest-mass energy eV / {m-nC?) 
and — eV/(m,rC^) respectively and the proton gravitational potential in units of the pion mass mp{e'^'^ — l)/m„. 
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FIG. 3: Left panel: electron number density for r > Rp normalized to its value at r = Rp. Right panel: electric field for r > Rp 
normalized to its value at r = Rp. We have also shown the behavior of the solution for two different eigenvalues of the general 
relativistic Thomas- Fermi equation (|lip close to the one which gives the globally neutral configuration (see 6] for details). 



(P/^)'^/(37r^), where Pf denotes the Fermi momentum of 
the j-specie. 

From Eq. ([28]) follows the Gibbs-Duhem relation 
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which can be rewritten as 
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Introducing the metric given by Eq. ([3]) the Einstein- 



Maxwell system of equations is 
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Using the Gibbs-Duhem relation (I5U1) the energy- 
momentum conservation equation (j33p can be rewritten 
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FIG. 4: Electric field (multiplied by lO") in units of the 
critical field Ec = nisC /{eh) ~ 10^ Volt/cm in the region 
r < Rn both for the Newtonian and the General Relativistic 
configurations. The central density of both systems is p(0) = 
3.94/9nuc where pnuc ~ 2.7 x 10^* g cm"'' is the nuclear density. 
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The Tolman isothermal condition [15| (see also jll|) 
demands the constancy of the gravitationally red-shifted 
temperature 

h -du = , or Too = e^/^T = constant , 

(36) 
which can be used into Eq. (1551) to obtain 

Y^ n,die''^^H^) + e{np~ne)dV = 0. (37) 



We now introduce the generahzed chemical potentials, 
or Klein potentials, for electrons Eg, protons Ep and neu- 
trons En 



Ee = c'^/Ve - rriec" - e^ , 
Ep = c^' Up — nipC + eV , 
E„. = e"^/2 
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which in the zero temperature limit are the generalized 
Fermi energies for electrons Ee = E^ , neutrons En = E^ 
and protons Ep — E^ introduced in Sec. II (see Eq. ([4])). 
Using Eqs. dSH]), dSlf and gO]), Eq. dST]) becomes 



y^ TiidEi = , 



(41) 



which leads for independent and non-zero particle num- 
ber densities n^ 7^ to the constancy of the Klein poten- 
tials ([38l) -([40 |l for each particle-species, i.e. 



Ee = e'^/^/ie - meC^ - cV = constant , (42) 
Ep = e^/'^iip ~ nipC^ + eV ^ constant , (43) 
En — e'^/^/in — ninC^ — Constant . (44) 

In the zero temperature limit the constancy of the 
Klein potential of each particle-specie becomes the con- 
stancy of the generalized Fermi energies introduced in 
Sec. [ni (see Eqs. (O and ©). This is a crucial point 
because, as discussed in [6|, the constancy of the gener- 
alized Fermi energies proves the impossibility of having a 
self-consistent configuration fulfilling the condition of lo- 
cal charge neutrality and /3-equilibrium (see e.g. Fig. [1}. 
Further, as shown in 16], the constancy of the Klein po- 
tentials holds in the more general case when the strong 
interactions between nucleons are taken into account. 

Therefore, introducing the new dimensionlcss variables 
Vi — P-i/i^T) and /3.i — kT/{mic'^), the new set of 
Einstein-Maxwell-Thomas-Fermi equations generalizing 
the system ([6|)-- (fTT|) to the case of finite temperatures 
is 
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where Eq. (j49p is the condition of /3-equilibrium between 
neutrons, protons and electrons, and the number density 
of the i-specie is given by 

"^ = ^-^^f^'^'(Pl/2 + ft^3/2) , (52) 

where we have introduced the relativistic Fermi-Dirac in- 
tegrals of order j 

F; = FAv.,^d-f"-^^j^j^d.. (53) 

We have integrated numerically the system of equa- 
tions (H5|) - (|5T|) for given temperatures Too ^ 0. As ex- 
pected, the results are both qualitatively and quantita- 
tively similar to the ones obtained with the degenerate 
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FIG. 5: Electron number density for r > Rp normalized to its 
value at r = _Rp both for T = K (degenerate case) and for 
a finite temperature of Too = 2.3 x 10^ K. 



approximation. The largest difference we found is at the 
surface boundary of the configuration, where due to the 
low density finite temperature effects are more effective. 
In Fig. [5] we compare the electron density for r > Rp 
in the degenerate and in the non-degenerate case for 
Too — 2.3 X 10^ K. For distances r < Rp the results are es- 
sentially the same as in the degenerate case. In the region 
r << Rn at large densities > pnuc — 2.7 x 10^* g/cm^, 
the electrodynamical properties of the configuration i.e. 
Coulomb potential and electric field remain unperturbed 
even for very large temperatures Too '^ 10^^ K. This is 
due to the fact that thermal effects are largely compen- 
sated by the gravitational potential as given by Eq. (|55| ; 
the Coulomb interaction is not involved in this balance 
and is not affected by the thermal energy. 



non-relativistic c —> oo limit of the general relativis- 
tic Thomas-Fermi and Einstein-Maxwell equations ([6])- 
([TTj) . The numerical integration of the Newtonian equa- 
tions shows that the gravito-electrodynamic structure ev- 
idenced in JI (see also Sec.|ll]) is not of general relativistic 
nature but is already present in the Newtonian regime. 
However, in view of the large quantitative discrepancies 
of the Newtonian regime in the description of other neu- 
tron star properties like mass and radius (see e.g. Fig. 3 
in 17]), a general relativistic treatment is mandatory in 
any astrophysical consideration. 

2) It has been also presented the extension of the pre- 
vious treatment [6( to finite temperatures. Although the 
thermal energy expected to be stored in old neutron stars 
with surface temperatures ^-^ 10^ K [ij,[ij is much larger 
than the internal Coulomb energy (see Sec. lIVp . still the 
electromagnetic structure (see Fig. H]) is unaffected by 
the presence of the thermal component. Physically this 
is due to the very large Fermi energy of the neutrons ^ 1 
GeV, of the protons ~ 10 MeV and of the electrons ^ 0.1 
GeV, as can be seen from Eq. ([23|) . 

3) More generally, we have given the explicit demon- 
stration of the constancy throughout the configuration of 
the Klein potentials of each species present in the system 
in the more general case of finite temperatures. This gen- 
eralizes the condition of the constancy of the general rel- 
ativistic Fermi energies derived in the special case T = 
in [6]. 

All these results are relevant to the generalization of 
the Feynman-Metropolis- Teller treatment of compressed 
atoms to relativistic regimes in presence of thermal ef- 
fects and they will necessarily also apply in the case of 
additional particle species and of the inclusion of strong 
interactions in neutron stars llal. 
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V. CONCLUSIONS 

In this article we have addressed three new aspects of 
the description of a self gravitating system of neutrons, 
protons and electrons in /3-equilibrium: 

1) We have presented the Newtonian limit of the treat- 
ment by taking the weak field approximation and the 
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